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ABSTRACT 



To improve the planet detection efficiency, current planetary microlensing 
experiments are focused on high-magnification events searching for planetary 
signals near the peak of lensing light curves. However, it is known that cen- 
tral perturbations can also be produced by binary companions and thus it is 
important to distinguish planetary signals from those induced by binary com- 
panions. In this paper, we analyze the light curves of microlensing events 
OGLE-2007-BLG-137/MOA-2007-BLG-091, OGLE-2007-BLG-355/MOA-2007- 
BLG-278, and MOA-2007-BLG-199/OGLE-2007-BLG-419, for all of which ex- 
hibit short-term perturbations near the peaks of the light curves. From detailed 
modeling of the light curves, we find that the perturbations of the events are 
caused by binary companions rather than planets. From close examination of 
the light curves combined with the underlying physical geometry of the lens sys- 
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tem obtained from modeling, we find that the short time-scale caustic-crossing 
feature occurring at a low or a moderate base magnification with an additional 
secondary perturbation is a typical feature of binary-lens events and thus can be 
used for the discrimination between the binary and planetary interpretations. 

Subject headings: gravitational lensing - planetary systems 



1. Introduction 



The microlensing signal of a planet is characterized by a short-term perturbation to the 
smooth standard single- l ens light curve of the primary-induced lensin g event occurring on a 
background source star jMao & Paczvhskilll99ll ; bould & Loeblll992h . The duration of the 
perturbation is several days for a gas giant and several hours for an Earth- mass planet, while 
the typical cadence of lensing surveys is roughly a day As a result, it is difficult to detect 
planets from the survey observations alone. To achieve the cadence required to detect short- 
term planetary signals, current planetary lensing searches are being co nducted by combining 



survey and follow-up observat ions, in which the survey observations (lOGLE: Udalski et al. 



20031 ; IMOA: Bond et al.ll2002l ) aim to maximize the number of detections of lensing events 



by monitoring a large area of sk y and follow-up observations (IPLANET: Albrow et al.lll998 



Micro-FUN: Ghosh et al.l 120041 ) are focused on intensive monitoring of the events detected 
by the survey observations. However, the number of telescopes available for follow-up ob- 
servations is far smaller than would be needed to monitor all the detected events from the 
survey observations. As a result, follow-up observations are selectively conducted for events 
that can maximize the planet detection probability. Currently, the highest priority is given 
to high-magnification events. For these events, the planet detection efficiency is high because 
the source trajectory always passes close to the region of central perturbations. In addition, 
the perturbation occurs near the peak of the light curve and thus follow-up observations can 
be prepared in advance. 

Although the strategy to detect central perturbations has an important advantage of 
high sensitivity to planets, it also has a shortcoming of difficulty in the interpretation of 
observed signals. One important cause of this difficulty is that the perturbation near the 
peak of a lensing light curve can be produced not only by a planet but also by a wide or a 
close binary companion with a mass roughly equal to that of the primary. Fortunately, the 
perturbat ion patterns produced by the p lanetary and the binary companion are intrinsically 
different ( IHan fc Gaudill2008l ; lHanll2009l ). and thus it is in principle possible to discriminate 
between the planetary and binary interpretations. However, this discrimination usually 
requires detailed modeling of the light curve, which demands a time-consuming search for a 
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solution of the lensing parameters in the vast space of many parameters. Therefore, simple 
diagnostics that can be used for the discrimination between the two interpretations will be 
useful not only for the immediate identification of the nature of the perturbation but also for 
the preparation of observational strategies for better characterization of planetary systems. 

In this paper, we present the results of the analysis of three binary-lens events with very 
strong short-term signals near the peaks of the light curves. From the investigation of the 
light curves combined with the underlying physical geometry of the lens system, we find and 
present several features in lensing light curves that can be used as diagnostic tools for the 
discrimination between the planetary and binary interpretations of observed events. 

The paper is organized as follows. In section 2, we briefly describe the general features of 
central perturbations produced by a planet and a binary companion. In section 3, we discuss 
the data used in our analysis. In section 4, we explain the modeling procedure and present 
the best-fit lensing parameters of the individual lensing events obtained from the modeling. 
In section 5, we present the features in the light curves that characterize the binary nature 
of the events and discuss the usability of the features in distinguishing the planetary and 
binary interpretations. We summarize the results and conclude in section 6. 

2. Planetary and Binary Central Caustics 

Lens systems composed of multiple masses can produce strong perturbations when the 
trajectory of a source star passes close to a caustic. The caustic represents the source star 
position at which the magnification of a point source becomes infinite. For lenses consisting 
of two masses, the set of caustics form closed curves each of which is composed of concave 
curves that meet at cusps. 

Binary lenses can have one, two, or three sets of caustic depending on the separation 
between the lens components. If the two masses are separated by approximately an Einstein 
radius, then there is a single six-cusp caustic. If the separation is substantially closer than 
the Einstein radius (close binary), there is a single four-cusp central caustic and two small 
outlying three-cusp caustics. If the masses are separated by substantially more than the 
Einstein radius (wide binary), there are two four-cusp caustics each of which is associated 
with a member of the binary. The short-term central perturbation of a binary-lens event can 
occur in two different ways. The first case is when the trajectory of a source star passes close 
to the central caustic of a close binary. The other case happens when the trajectory passes 
close to one of the two caustics of a wide-separation binary. For both wide and close binaries, 
the central caustics have a diamond shape although there is some variation depending on 
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the separation and the mass ratio between the two binary components. 



For planetary lensing, the central caustic refers to a small four-cusp caustic located close 
to the primary lens. The size of the planetary central caustic as measured by the separation 
between the tw o cusps located on the primary-planet axis is related to the planetary lensing 
parameters by (IChung et al.ll2005l ) 



iq 



, s — s 



-1\2 ' 



where s represents the projected primary-planet separation normalized by the Einstein radius 
corresponding to the total mass of the lens system, #e, and q is the planet /primary mass 
ratio. While the size of the central caustic depends on both the mass ratio and the separation, 
its shape depends solely on the separation. When the shape is quantified as the ratio of the 
vertical width Arj r , to the horizon tal width A£ c , the width ratio is related to the planetary 
separation by (IChung et al.ll2005l ) 



Ar/ C 



sin 



(2) 

A£ c (s + s- 1 -2cos0) 2 ' V ; 

where cos0 = (3/4) (s + — [(32/9)(s + s -1 ) 2 ] 1 / 2 }. For the range of the planetary 

separations within which the planet detection efficiency is not ne gligible, th e width ratio of 
the planetary central caustic is substantially smaller than unity (jHanll2009l ). implying that 
the caustic is elongated along the primary-planet axis. The caustic has an arrow-head shape 
where three cusps are located close to the primary and the other cusp is located at the tip 
of the arrow pointing toward the direction of the planet. For a given mass ratio, a pair of 
central caustics with sep arations s and s^ 1 are identical to the first order of approximation 
jDominiklll999l : [Anlhoosh . 



3. Observations 

The events we analyze include OGLE-2007-BLG-137/MOA-2007-BLG-091, OGLE-2007- 
BLG-355/MOA-2007-BLG-278, and MOA-2007-BLG-199/OGLE-2007-BLG-419. All these 
events were observed in the 2007 microlensing observation season. The light curves of the 
individual events are presented in Figure [1] - Figure [3j They show a common feature of 
strong perturbations occurring near the peaks of the light curves. Due to the location of the 
perturbations similar to those of central perturbations induced by planets, all events were 
alerted by survey observations and intensively monitored by follow-up observations. 

The lensing event OGLE-2007-BLG-137/MOA-2007-BLG-091 was alerted on 2007 April 
15 by the OGLE collaboration as a possible event with an anomaly near the peak, using the 
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1.3 m Warsaw telescope of Las Campanas Observatory in Chile. In response to the alert, 
follow-up observations were conducted by the /iFUN team with various telescopes including 
the 1.3 m of CTIO in Chile, the 0.4 m of Auckland Observatory, the 0.4 m of Vintage Lane 
Observatory in New Zealand, and the 1.0 m of Mt. Lemmon Observatory in Arizona, United 
States. The MOA collaboration independently discovered the event by using the 1.8 m of 
Mt. John in New Zealand. From these observations, a total of 457, 1119, 14, 65, 2, and 2 
images were obtained from Las Campanas, Mt. John, CTIO, Auckland, Vintage Lane, and 
Mt. Lemmon, respectively. 

The event OGLE-2007-BLG-355/MOA-2007-BLG-278 was detected independently by 
the OGLE and MOA survey grou ps and its anom alous behavior was noticed by the SIG- 



NALMEN, a software developed by iDominikl (120071) for ea rly detection of anomalies and now 



integrated part of the ARTEMiS system (lDominikll2008l ). on 2007 July 20. Soon after, the 



//FUN and PLANET/RoboNet teams conducted intensive follow-up observations by using 
the 1.3 m of CTIO, the 0.4 m of Auckland Observatory, the 0.36 m of Farm Cove Obser- 
vatory in New Zealand, the 0.28 m of Southern Stars Observatory in Tahiti, the 0.36 m of 
Bronberg Observatory in South Africa, 1.54 m Danish Telescope of La Silla Observatory in 
Chile, the 1.0m of SAAO in South Africa, 2.0 m Faulkes S. (FTS) Telescope in Australia, 2.0 
m Faulkes N. (FTN) Telescope in Hawaii, and 2.0 m Liverpool Telescope (LT) in La Palma, 
Spain. Due to the prompt response to the alert, the perturbation was very densely resolved 
by the follow-up observations despite the fact that the perturbation lasted only about 1 day. 
The data consist of 591, 1325, 187, 21, 39, 3, 43, 40, 34, 3, 27, and 7 images obtained from 
Las Campanas, Mt. John, CTIO, Auckland, Farm Cove, Southern Star, Bronberg, La Silla, 
SAAO, FTS, FTN, and LT, respectively. 

The event MOA-2007-BLG-199/OGLE-2007-BLG-419 is a long time-scale event. The 
PLANET/RoboNet team first announced the anomalous behavior of the event on 2007 July 
28 and conducted follow-up observations by using the 1.54 m Danish Telescope, the 1.0 m of 
Canopus Observatory, the 0.6 m of Perth Observatory in Australia, and the 1.0m of SAAO. 
The /xFUN collaboration also observed the event by using the 1.3m of CTIO. The event 
shows a clear signature of caustic crossing during HJD ~ 2454316 - 2454322. The data 
consist of 1233, 4108, 126, 56, 52, 130, and 15 images obtained from Las Campanas, Mt. 
John University, La Silla, Canopus, Perth, SAAO, and CTIO, respectively. 



For the analysis of the light curves, we use the photometry results that were reduced by 
the individual groups using their own photometry codes. 
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4. Modeling 

Binary-lens modeling requires inclusion of various parameters. The first set of three 
parameters is related to the geometry of a single- lens event. These include the Einstein time 
scale, the time of maximum magnification, t , and the lens-source separation normalized 
by the Einstein radius 6*e at the time of maximum magnification, uq. For the description 
of the perturbation, an additional set of binary parameters is required. These parameters 
include the mass ratio q and separation s between the lens components and the angle of 
the source trajectory with respect to the binary axis, a. In addition, the normalized source 
radius, p* = 9+/9 E , is needed to describe the lensing magnification whenever the angular 
radius of the source star 6+ plays an important role, such as the moment of a caustic crossing. 
For some events with long time scales, it is required to include the parallax parameters 7Te,n 
and 7Te,e to account for the deviation of the source motion from r ectilinear motion caused by 



the parallactic motion of the Earth around the Sun (lGouldll200d ). The microlens parallax is 
defined by n E = AU/f E , i.e., the ratio of the radius of the Earth's orbit around the Sun to 
the Einstein radius projected on the observer plane, r E . The subscripts "iV" and U E" refer 
to the components of the parallax vector projected on the sky in the north and east celestial 
coordinates. 

To determine the model parameters, we conduct a massive search for solutions over 
a broad range of mass ratios and separations. For this, we hold the binary parameters s, 
q, and a fixed at a grid of values, while the remaining parameters (£e, uo, to, and p*) are 
allowed to vary so that the model light curve results in minimum x 2 a ^ each grid point. We 
use a Markov Chain Monte Carlo method for the \ 2 minimization. Once the x 2 minima of 
the individual grid points are determined, the best-fit model is obtained by comparing the 
X 2 values. Although a brute search in the space of the binary parameters requires a large 
computation time, it is needed to investigate the existence of possible degenerate solutions. 

For the computation of lensi ng magnification includi ng; the finite-source effect, we use the 



magnification map" technique (IDong et al.ll2006l . 120091 ). In this technique, a magnification 



map is constructed for a given set of s and q by using the ray-shooting method. Based on 
this map, the lensing magnification corresponding to a location of a source star on the map is 
computed by comparing the number density of rays within the source radius with the density 
on the lens plane. On the basis of this scheme, the computation is accelerated in several 
ways. First, the algorithm is designed to save computation time by limiting the range of ray 
shooting on the image (lens) plane to a narrow annulus encircling the Einstein ring of the 
primary lens. This is based on the fact that only the rays from this region arrive at the central 
perturbation region on the source plane. Second, the algorithm keeps the information of the 
positions of all the rays arriving at the target in the buffer memory space of a computer so 
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that the information can be promptly used for fast computation of the magnification. Third 
finite-source magnification calculation based on ray-shoot ing is used only in the regi o n near 



the c austic and a simple hexadecapole approximation fjPejcha &: Heyrovskyi 120091 ; iGould 



20081 ) is used in other part of light curves. 



In Figure [T] - Figure [3j we present the model light curves superposed on the data of 
the individual events. The geometry of the source trajectory with respect to the caustic 
location for each event is presented in the inset of each figure. The coordinates of the inset 
are centered at the center of magnification around which the single-lens magnification best 
describes that of the binary lensing. The abscissa and ordinate are parallel with and normal 
to the binary axis, respectively. In Table ??, we also present the values of the individual 
lensing parameters^ We note that modeling the light curve of the event MOA-2007-BLG- 
199/OGLE-2007-BLG-419 with the basic binary parameters does not yield a satisfactory 
model and requires inclusion of additional parallax parameters. As a result, the trajectory 
of the event is curved. It is found that the determined mass ratios of the three lenses range 
from 0.24 to 1.0, implying that the perturbations are produced by binary companions not 
planets. It is also found that the short strong perturbations are caused by the crossing the 
tip of the caustic. 

We note that including the parameters of higher-order effect, p* and 7r E) is necessary 
to obtain a good fit to the data. But the data reductions, which comes directly from the 
pipelines, might not be adequate to measure these subtle parameters accurately. Since these 
parameters are irrelevant to the main scientific purpose of the paper, and since refined data 
reductions are typically very time-consuming, we present these higher-order parameters only 
for completeness. 



5. Binary Features 

The light curves of the events we analyzed are mainly characterized by the perturba- 
tion near the peak similar to the central perturbation caused by a planet. However, close 
examination of the perturbations reveals several features that convey the binary nature of 
the perturbation. In this section, we discuss these features. 

The first feature that gives a clue as to the nature of the events is the strength of the 



x For the event MOA-2007-BLG-199/OGLE-2007-BLG-419, we present two sets of solutions resulting from 
the known "+wo — u o" degeneracy in parallax determination ( Smith et alj|2003l) . We find tha t there exists 
no "jerk-parallax " degeneracy, which is another known degeneracy in parallax determination (jGouldl 12004 ; 



Park et al.ll2004l) . 
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main perturbation. For all three events, the deviations from the smooth base of the light 
curve exceed 3 magnitudes. This strength suggests that the perturbations are caused by 
caustic crossings. 

The second feature is the duration of the perturbation. For all events, the duration of 
the main perturbation is very short compared to their time scales. The short duration of the 
perturbation can be explained either by the caustic being tiny or by the source trajectory 
crossing the tip of a cuspy part of a caustic with a non-negligible size. 

The third feature is the base magnification just before and after the perturbation. The 
base magnification of a caustic-crossing event is important for the characterization of the 
lens system because it gives a strong constraint on the size of the caustic. For the event 
OGLE-2007-BLG-355/MOA-2007-BLG-278, the base magnification is only A ~ 2 (corre- 
sponding to the normalized lens-source separation of u ~ 0.5). This implies that the caustic 
size is half as large as the Einstein radius and thus the binary origin of the perturbation 
can be immediately noticed. The events OGLE-2007-BLG-137/MOA-2007-BLG-091 and 
MOA-2007-BLG-199/OGLE-2007-BLG-419 have moderate base magnifications of A ~ 10 
(corresponding to u ~ 0.1). Considering that the size of the central caustic produced by 
a giant planet with a mass ratio q ~ 10~ 3 located around the Einstein ring of the primary 
with a separation s ~ 1.3 is A£ c ~ 0.01, the chance that the main perturbations of these 
events are caused by planets is not high. However, one still cannot rule out the possibility of 
a big central caustic produced by a planet located very close to the Einstein ring (resonant 
planet). 

The last feature is the additional secondary perturbation. The central caustic of a 
resonant planet can become large and thus planetary perturbations can occur at a moderate 
base magnification. Figure H] shows the magnification maps of a resonant planetary system 
in the region around the central caustic. To produce a planetary perturbation at a low 
magnification, the source trajectory should pass close to the cusp corresponding to the tip 
of the arrow-head shaped planetary central caustic such as the one shown in Figure HI In 
this case, it is difficult to produce an additional secondary perturbation because all the other 
cusps are remotely located on the opposite side of the caustic close to the primary. On the 
other hand, the central caustic produced by a binary companion has a diamond shape and 
thus the geometric chance to produce a secondary perturbation by passing an adjacent cusp 
of the caustic is high as demonstrated in Figure HI Therefore, the existence of secondary 
perturbation strongly advocates the binary interpretation of the events. 

In summary, central perturbations that are characterized by short time-scale caustic- 
crossing features occurring at a low or a moderate base magnification with an additional 
secondary perturbation are caused by binary companions. Considering that three events 
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were detected during a single observation season, events with such a perturbation feature 
are common and the feature can be used as a diagnostic for the discrimination between the 
binary and planetary interpretations of perturbations. 

6. Conclusion 

We analyzed the light curves of three microlensing events with short-term perturbations 
near the peaks of the light curves. From detailed modeling of the light curves, we found 
that the perturbations of the events are caused by binary companions rather than planets. 
From close examination of the light curves combined with the underlying physical geometry 
of the lens system obtained from modeling, we found that the short time-scale caustic- 
crossing feature occurring at a low or a moderate base magnification with an additional 
secondary perturbation is a typical feature of binary-lens events and thus can be used for 
the discrimination between the binary and planetary interpretations. 
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Fig. 1.— Data and best-fit model of the microlensing event OGLE-2007-BLG-137/MOA- 
2007-BLG-091. The left inset shows a blowup of the region of major perturbation. The 
right inset shows the geometry of the source trajectory with respect to the caustic location. 
The coordinates are centered at the 'center of magnification' of the lens system (see the 
definition in the text) and the abscissa and the ordinate are parallel with and normal to 
the binary axis, respectively. Heavier mass is on the left. The temperature scale represents 
the magnification where brighter tone represents higher magnification. The bottom panel 
represents the residual from the model light curve. 
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Fig. 2.— Data and best-fit model of OGLE-2007-BLG-355/MOA-2007-BLG-278. Notations 
are same as in Fig. [TJ 
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Fig. 3.— Data and best-fit model of MOA-2007-BLG-199/OGLE-2007-BLG-419. Notations 
are same as in Fig. [1] except that there are two insets of source trajectory geometry for the 
two degenerate models. The presented model light curve is for "+Wo" solution (see Table 
??), which has smaller Ay 2 by ~ 90. 
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Fig. 4. — The magnification maps of a resonant planetary system (upper left panel) and 
a close binary system (upper right panel) in the regions around the central caustics of 
the individual systems. The lower panel shows the light curves resulting from the source 
trajectories marked in the upper panels. 



